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NOMENCLATURE 


a blade lift-curve slope 

ao blade coning angle measured from hub plane, rad (or deg) 

a\ longitudinal first-harmonic flapping coefficient measured from hub plane 
and in wind-hub system, rad (or deg) 

a^g longitudinal first-harmonic flapping coefficient measured from hub plane 
and in hub-body system, rad (or deg) 

lateral cyclic pitch measured from hub plane and in wind-hub system, 
rad (or deg) 

A^g lateral cyclic pitch measured from hub plane and in hub-body system, 
rad (or deg) 

b| lateral first-harmonic flapping coefficient measured from hub plane and 
in wind-hub system, rad (or deg) 

b^g lateral first-harmonic flapping coefficient measured from hub plane and 
in hub-body system, rad (or deg) 

longitudinal cyclic pitch measured from hub plane and in wind-hub 
system, rad (or deg) 

B|g longitudinal cyclic pitch measured from hub plane and in hub-body system, 
rad (or deg) 

c blade chord, m 


B 

K 


B 


B 


w 


flapping hinge offset, m 

blade moment of inertia about flapping hinge, kg-m^ 

pitch-flap coupling ratio, = tan 63 

flapping hinge restraint, N-m/rad 

blade mass moment about the flapping hinge, kg-m 

number of blade 

aircraft roll rate in hub-body system, rad/sec 
aircraft roll acceleration, rad/sec^ 

aircraft roll rate in wind-hub system, rad/sec, = p cos 6^ + q sin 6^ 


iii 


ratio of flapping frequency to rotor system angular velocity 
aircraft pitch rate in hub-body system, rad/sec 


aircraft pitch acceleration, rad/sec^ 

aircraft pitch rate in wind-hub system, rad/sec, 
= -p sin + q cos 6^ 


w 


w 


radial station of the blade element measured from the flapping 
hinge, m 

rotor radius, m 

uniform induced velocity, m/sec 
true airspeed, m/sec 


nondimensional radial station of the blade element, x ^ 


e + r' 


wind-hub system 


hub-body system 


hub plane angle of attack, deg 

blade flapping angle measured from hub plane, rad (or deg) 


A dR 

derivative of flapping with respect to time, g £ -~ 

dt 

second derivative of flapping with respect to time. 


dt^ 


differential collective flapping (only for even-bladed rotors) 
nth order longitudinal cyclic flapping 
nth order lateral cyclic flapping 
collective flapping (coning) 


multiblade flapping coordinates 
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w 


Y 


9 




0 


t 


A 


y 

p 


c 


Q 


rotor sideslip angle, that is, the angle between and x^, deg 
Lock number, = — 

lateral control displacement, cm 

longitudinal control displacement, cm 

e 

R 

blade pitch angle measured from hub plane, 

0 = 0Q - cos ip - sin jp + K0^ - Kj3? f^d (or deg) 

damping ratio 

blade-root collective pitch measured from hub plane, rad (or deg) 

total blade twist (tip with respect to root), deg 

V sin a - vj 
inflow ratio, 

QR 

V r*n c: rv 

advance ratio, — 

0 u R. 

air density, kg/m^ 
rotor solidity ratio 
control advance angle, deg 

azimuth angle measured from downwind in the sense of rotor rotation, 
rad (or deg) 

azimuth angle measured from -x^ in the sense of rotor rotation, 
rad (or deg) 

undamped natural frequency, rad/sec 
rotor system angular velocity, rad/sec 


Subscript 

s.s. steady state 
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EFFECTS OF PRIMARY ROTOR PARAMETERS ON FLAPPING DYNAMICS 


Robert T. N. Chen 
Ames Research Center 


SUMMARY 


An analysis has been made to study the effects of flapping dynamics of 
four main rotor design features that influence the agility, stability, and 
operational safety of helicopters. The parameters include flapping hinge 
offset, flapping hinge restraint, pitch-flap coupling, and blade Lock number. 
First, the flapping equations of motion are derived that explicitly contain 
the design parameters. The dynamic equations are then developed for the tip- 
path plane, and the influence of individual and combined variations in the 
design parameters determined. 

The analysis includes a study of the steady-state flapping response with 
respect to control input and aircraft angular rate which leads to a feedforward 
control law for control decoupling through cross feed, and a feedback control 
law to decouple the steady-state flapping response. The condition for achiev- 
ing perfect decoupling of the flapping response due to aircraft pitch and roll 
rates without using feedback control is also found for the hover case. 

The analysis also indicates that the frequency of the regressing flapping 
mode of the rotor system can, for some designs, become low enough to require 
consideration in the assessment of handling characteristics (less than 2 Hz). 
This occurs for rotors with a large effective hinge offset and with heavy rotor 
blades . 


INTRODUCTION 


Continual improvement is sought in the inherent agility and stability of 
the helicopter to expand its operational capabilities while providing adequate 
safety margins. Recently, considerable emphasis has been placed on helicopter 
main rotor design features to enhance control moments and stability. Improve- 
ments in these areas are required, for instance, to permit single pilot commer- 
cial operation under instrument flight rules (IFR) , and to provide the agility 
and maneuverability essential for nap-of-the-earth flight. The important 
design features include several that increase control moments by providing a 
direct hub moment in addition to that generated by thrust vector tilt. These 
features are the hingeless rotor (ref. 1) and the stiffened flapping hinge for 
the teetering rotor system (ref. 2). To improve the stability in forward 
flight, particularly with the hingeless rotor, consideration has been given to 
the use of mechanical means such as pitch-flap coupling (ref. 3). Safety of 
flight considerations, particularly for single-engine helicopters, leads to the 
high rotor inertia concept (ref. 4) to eliminate or reduce the dead man’s 
height -velocity region. 


Specific direct requirements may be achieved through the use of these 
main rotor design features. However, side effects, such as inter-axis cou- 
pling (ref. 5) and control response lag due to rotor dynamics, can adversely 
affect handling qualities. Insight on the potential impact of these design 
features on handling characteristics can be gained through a detailed examina- 
tion of their effects on flapping dynamics. The analytical study reported 
herein provides this insight through consideration of the effects of large 
variation of flapping hinge offset, flapping hinge restraint, blade Lock 
number, and pitch-flap coupling on the flapping dynamics. Special emphasis 
was placed on determining the effects of these four design parameters on the 
tip-path plane dynamics and on the coupling of the longitudinal and lateral 
flapping due to control input and due to pitch and roll motion. 


FLAPPING EQUATION OF MOTION 


The flapping equation of motion was derived explicitly for a rotor system 
with the four parameters of interest, as shown in figure 1: the stiffness of 

the flapping hinge, K^; the effective hinge offset, e; the blade Lock number, 
y; and the pitch-flap coupling, K^. 

To develop analytical expressions that may provide a better insight into 
the parametric effects, many simplifications and assumptions similar to those 
used for the ’’classical'* equations (ref. 6) were used. They are: 

1. Rotor blade was rigid in bending and torsion, and the twist of the 
blade was linear. 

2. Both the flapping angle and inflow angle were assumed to be small and 
the analysis utilized a simple strip theory. 

3. The effects of the aircraft motion on the blade flapping were limited 
to those due to the angular acceleration p and q, the angular rate p, q, and 
the normal acceleration. 

4. The reversed flow region was ignored and the compressibility and stall 
effects disregarded , 

5. The inflow was assumed to be uniform and no inflow dynamics were used, 

6. The tip loss factor was assumed to be 1. 

Because of these assumptions and simplifications, the results of the anal- 
ysis were valid only for a limited range of flight conditions. Nevertheless, 
a previous study (ref. 7) has shown that this type of analysis is valid for 
stability and control investigations of the rotorcraft up to an advance ratio 
of approximately 0.3. 

The development of the flapping equation of motion as described in appen- 
dix A was rather elementary and straightforward. The wind-hub coordinate 
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system, Xg, yg, Zg, as depicted in figure 2, was employed as the basic frame 
of reference. In this coordinate system, Xg, yg lie on the hub plane of the 
rotor, and the Xg-Zg plane contains the relative wind. 


For nonteetering rotor systems, the flapping equation for a single blade 
was obtained by summing the moments due to aerodynamic forces acting on the 
blade, the centrifugal force, Coriolis, inertia, restraint forces, and the 
blade weight about the flapping hinge. The terms containing the hinge offset 
were retained up to the second power. The result as derived in appendix A is 
shown in equation (1). 
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where 


P-" = 1 + 




‘6 


I 

V" 


8 V 3 V 


i = 1, 2, 3, . . ., N 
= p cos 6^ + q sin 

" "P q ^w 

For a two-bladed teetering rotor, the flapping equation was derived from equa- 
tion (1) using the following constraints: 


e = 0 

The result is shown in the following: 
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where 


K 


= 1 + 




3 , Ki^ 


To gain a better insight into the dynamics of the flapping motion, equa- 
tion (1) will be transformed into a nonrotating coordinate system using the 
multibladed coordinate transformation (ref. 8), that is. 


n=i 


3. = 3o + 3, (-1)^ P n^j . + 3 sin nil^ . 

1 cl nc 1 ns 1 

i^lj 2, 3, . , »,N 


(3) 


where 
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be denoted, respectively, by 


N odd 
N even 


(B, = 0, if N odd) 
cl 


_B - ( 3o j 3ie > 3is » 62 c 5 ^ 2 s > • • ■’ ^ ( 1 / 2 ) (N- 1 ) c ’ ^ ( 1 / 2 ) (N- 1 ) s^ 

T 

k (6o.3ic»Sis> • • •» (N-2)c’ 1 / 2 ) (N- 2 )s’ 


for N = odd 
for N = even 


Then, as shown in appendix B, the flapping equation (1) in the rotating coor- 
dinate system may be transformed into the nonrotating coordinate system to 
yield : 


6 + DB + K 6 = f 


(A) 


Tables 1 and 2 show the flapping equation in the nonrotating coordinate 
system for the three-bladed and four-hladed rotor systems, respectively. From 
these tables, the following observations may be made: 

(a) In nonrotating coordinates, the flapping equations also contain peri- 
odic coefficients in forward flight. The basic frequency of the periodicity 
is directly related to the number of blades in the rotor; the basic frequency 
is 3 per rev for the three-bladed rotor and 2 per rev for the four-bladed 
rotor. The highest frequency in the periodic terms is N per rev for an 
N-bladed rotor. In general, it can be shown (as previously stated in ref. 9) 
that for 3, 5, 7, . . . bladed rotors, the basic frequency is 3, 5, 7, . . ., 
respectively; for 4, 6, 8, . . . bladed rotors, the corresponding basic fre- 

quency is 2, 3, 4, . . .. 
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(b) The amplitudes of the periodic terms in the damping matrix are func- 
tions of Y» £. and y; they are independent of K and . At a given advance 
ratio, the lower the Lock number of the blade and the larger the hinge offset, 
the smaller the amolltudes of the periodic terms in the damping matrix. 

(c) The maximum magnitude of the periodic coefficients in the stiffness 
matrix are functions of y, e, and . Again the hinge restraint has no direct 
impact on the periodic terms. At a given advance ratio, a decrease in the 

hock number will decrease the amplitude of the periodic terms in the 
stiffness matrix. However, the combined effect of e and is more compli- 
cated. For =0, an increase in the hinge offset will reduce the effect of 

the periodic terms. 

. parametric effect on the periodic terms in the forcing functions 

is similar to that on the damping terms for a given set of control positions, 
®o» ^ic* ®ic» the twist of the blade 0^. 

The above observations may provide a better insight in assessing the qual- 
ity of approximation when the time varying flapping equations (expressed in 
the nonrotating coordinate system) in forward flight are simplified to time 
invariant system of equations by dropping the periodic terms. At y = 0, of 
course, there will be no periodic terms in these equations. 

With the periodic terms dropped, the first three equations in table 2 for 
N - 4 collapse to those in table 1 for N = 3; the differential coning equa- 
txon (the fourth equation) becomes uncoupled and with zero forcing function. 

The set of the first three equations is identical to the first-harmonic approx- 
imation also known as the "tip-path plane" equation as discussed in the next 
section. 


TIP-PATH PLANE DYNAMICS 


The mathematical procedure for approximating the flapping equation by per- 
forming the multiblade coordinate transformation and then dropping the periodic 
terms as discussed in the previous section is identical to the classical 
method of approximating the flapping by the first-harmonic terms with time 
varying coefficients, that is, 

p(t) =ao(t) - aj^(t)cos tit-bi(t)sini|; (5) 

Equating, respectively, the constant term, and the terms with sin ij; and 
cos iij in the equation (1) using (5), yields the tip-path plane dynamic 
equations : 


^ + Da. + ica^ = f 

The damping matrix D, the stiffness matrix K, and the forcing function f 
are shown in table 3. This set of equations is identical to that in table 1 
with periodic terms dropped. The identity is readily seen by noting that 
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3o = 6o» ai = -3ic> and bi = “3is» performing a similarity transforma- 

tion for the dynamic equations. 

If the ”tip~path plane” representation procedure described above is 
extended to approximate the flapping equation for the teetering rotor, equa- 
tion (2) becomes 



and 

B(t) = ao - ai(t)cos - bi(t)sin 4' (8) 

Note that in equation (8), which approximates the flapping, the coning angle 
ao is a constant (a precone angle). The tip-path plane approximation for a 
two-bladed rotor is generally valid for only low frequency excitation; hence, 
the terms containing aircraft angular acceleration p and q in equation (7) 
may be omitted. 

Equation (7) is identical to the lower 2x2 block shown in table 3 with 
e set to zero. The quality of this approximation to the flapping equation of 
a teetering rotor system has not been fully evaluated. A preliminary numeri- 
cal assessment showed that for y > 0.4,^ the approximation in terms of sta- 
bility characteristics could deteriorate rapidly, especially when Ki = = 0 

^ further evaluation is therefore needed to fully explore its limitations. 


^Of course, at this condition, the validity of the model has already been 
violated, as described earlier in the report. 
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The dynamic characteristics of the tip-path plane as governed by the 
equations in table 3 will now be examined. There are three natural modes in 
the tip-path plane dynamics: coning, advancing, and regressing. Of these, 

the regressing flapping mode is most important concerning the effect of rotor 
dynamics on the handling characteristics of the rotorcraft. The regressing 
flapping mode is the lowest frequency mode of the three and it has a tendency 
to couple into the fuselage modes. The other two modes (the coning and the 
advancing) have higher undamped natural frequencies, respectively, on the order 
of rotating frequency and twice the rotating frequency of the rotor system; 
their impact on the rotorcraft handling characteristics is therefore consider- 
ably less significant. 


To gain insight into the effects of the parametric variations in the gen- 
eral rotor system on the modal characteristics, hover conditions will be con- 
sidered first. At hover, the coning equation is decoupled as evidenced in 
table 3. The coning mode has the undamped natural frequency and the 

damping ratio respectively, as follows: ^ 
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nc 
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OQ 
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eM^ 
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e + 2e^) 





(9) 


The undamped natural frequency and the damping ratio for the advancing and 
regressing modes are, respectively, given by 


nA 


= + 2P A - ) 


1/2 


nA 


In 


(10a) 


nR 


= + p2 - 2P/l - c;2y 


nR 


In 


/2 


(10b) 


Figures 3 through 5 show the effect of individual parametric variation in 
hinge restraint, and the pitch~flap feedback, on the eigenvalues 
of the advancing and regressing modes. The effect on the blade Lock number is 
also shown in these figures which illustrate that y is the most significant 
parameter in reducing the undamped natural frequency of the regressing flapping 
mode. If we concentrate on the region 5 0.5 (which corresponds to the 

region with undamped natural frequency of the regressing flapping mode of no 

2.5 Hz for a 300 rpm main rotor system) , the effect of a large hinge 
offset is also very significant in reducing up to y = 12. Thus, for a 
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small Y ^ large e, significant impact of the flapping dynamics on the 

handling characteristics of the rotorcraft can be expected. 

In forward flight, the coning equation is no longer decoupled and the 
modal characteristics of the tip-path plane dynamics become much more compli- 
cated. Figures 6 through 9 show the effect of advance ratio on the eigenvalues 
of the three modes for four sets of parameter values, respectively: 


K = e = Ki = 0 

P 

0 = 0.15 , Kp = Ki = 0 
— ^ = 0.225 , e = Ki = 0 

Ki = 0.225 , Kg = E = 0 


As Indicated in these figures, the tip-path plane dynamic modes are relatively 
insensitive to the variation in the advance ratio u (up to the validity of 
the math model, y = 0.3 to 0.4) for y < 8. With a large hinge offset, the 
eigenvalues of the TPP modes are almost invariant over the entire range of 
Lock number considered. 

As described earlier, the advancing mode of the TPP dynamics is on the 
order of 2 per rev (i.e., on the order of 10 Hz for a 300 rpm rotor system). 
Thus, the model will not be useful for real-time simulation if the computa- 
tional cycle time is larger than (or on the order of 0.050 sec for a 

300 rpm rotor system). If the budgeted cycle time cannot be reduced, some 
simplified models that will approximate the low frequency characteristics of 
the TPP dynamics may be desirable. A simple model that matches well with the 
equations in table 3 in the low frequency region may be obtained by simply 
setting a = 0, as suggested in reference 9. The result is 

a = - ( 11 ) 

Figures 10 through 12 show the comparison of characteristic roots of the 
reduced order model and the unreduced model at hover. For a forward flight at 
y =0.3, their characteristic roots are shown in table 4. The fidelity in the 
transient response is shown in figures 13 and 14. 


CONTROL PHASING AND FLAPPING RESPONSE DECOUPLING 


It is evident from equation (6) and table 3 that the flapping frequency 
will generally be different from the rotational frequency of the rotor system. 
Therefore, the maximiim flapping response to a cyclic-control input will no 
longer exhibit 90° lag in phase. Proper control phasing or control mixing will 
be required to achieve the desired flapping decoupling; a longitudinal control 
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input produces only a steady-state longitudinal flapping, and a lateral control 
input will produce only the lateral flapping response. 

To relate the required control phasing (or control mixing) to the rotor 
system parameters, consider first the steady-state solution of TPP variables 
at hover. From equation (6), the steady-state responses of the TPP variables 
ao> ai, and bi to the control inputs and the fuselage angular velocity p 
and q^ may be obtained using ^ 

^ (12) 


At hover, equation (12) becomes 



P = 1 + 


A- 


eM, 


+ 



4 

3 




In equations (13b) and (13c), the terms containing the fuselage angular accel- 
erations, p and q, have been dropped to be consistent with the steady-state 
solution of the TPP variables. From these equations, it is interesting to 
note that, for the special case P^ = 1, e = 0, the steady-state responses of 
the longitudinal and lateral flapping to the cyclic-control inputs reduce to 
the well known results: 


ai 


s . s. 


-B 


Ic 


b. = A 

s . s. Ic 
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Feedforward Control Law 


We will now develop the feedforward control law to achieve an appropriate 
control mixing for decoupling the steady-state flapping response to control 
displacements in the cockpit. To facilitate the flight- control system mech- 
anization, the development will be based on the hub-body coordinate system. 

In this coordinate system, the cyclic pitch will be denoted by and 

and the tip-path plane tilt by a^^ and b^g. They relate to those in the wind- 

hub system by the transformation: 


Ale 

= Ais 

cos 

^w 

— 
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sin 
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Bic 
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+ 
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Pw 
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sin 

Bw 


^is 

sin 



Bis 

cos 

Bw 


Consider the hovering flight first. At hover, the longitudinal and lat- 
eral flapping depend on both and but not on the collective pitch; 

the coning depends only on the collective pitch as evident from equations (13a), 
(13b), and (13c). Therefore, to decouple the steady-state flapping response, 
a control mixing of the following form will be satisfactory. 

“is , . “is . 

A L 

9B,g 3B,s 

®is " 96^ ^A 96^ \ 


In the above equations, 6^ and 5^ are, respectively, the lateral and longi- 
tudinal control displacement at the cockpit. The control gearings 
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96l 


for mixing the cyclic pitch may be chosen to satisfy a set of specified flap- 
ping response characteristics. Let the desired decoupled steady-state 
responses be 


where the parameters k] 
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( 14 ) 


Alternately, if the control gearings 9A^g/96^ and have been 

selected, decoupling the steady-state flapping may be achieved by appropriately 
phasing the cyclic-pitch control. The required control phasing (sometimes 
called control advance angle) in the sense of rotor rotation is given by 


= tan ^ [ (9B^g/96^) / ( 9 Ais/^^A^ ^ 


'[y[(1/8) - (e/3) + (eV4)]] 

Physically, the required control advance angle may be seen even more clearly 
by considering the flapping motion in the rotational coordinate. At hover, 
the flapping motion is governed by a second-order time invariant system with a 
natural frequency of VQ and a damping ratio of y/P[(l/16) - (c/3) H- (e^/8)], 
as evident from equation (1). At the excitation frequency Q, at which the 
cyclic pitch applies, the flapping response will have (90° - ())) lag in phase. 
For P = 1, when the flapping frequency is equal to the rotational frequency, 
then c}) = 0, as is well known. 


In forward flight, the coning is coupled to the cyclic flapping. To 
achieve a perfect decoupling a control law which utilizes a mixing of the 
cyclic stick to the collective pitch as well as the cyclic pitch is required. 
The development of the control law is given in appendix C. 

For a two-bladed teetering rotor system, a direct application of equa- 
tion (14) to equation (7) results in the following feedforward control law: 
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1 + (yK^yVie) 
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p2 - 1 + (3/16 )yK^P^' 

V ^ - 

1 - (yV2) 

Y 

1 + (3/2)y2 

K2 

1 + (3/2)y^_ 


(16) 


It is interesting to note that equation (16) reduces to equation (14) at hover. 

The decoupling control law equation (14) has been used previously in the 
nap-of- the-earth flight simulation (ref. 5) because of its simplicity in com- 
puter mechanization. The results have provided an approximation adequate for 
a speed range from hover up to an advance ratio of 0.25. 


Decoupling Flapping Response Due to Aircraft Angular Rate 

In the preceding discussion, attention has been focused on the decoupling 
of the flapping response to the control input. The other important aspect of 
decoupling the flapping response to the aircraft angular rate will now be dis- 
cussed, Unlike fixed-wing aircraft, for which pitch-roll coupling is rare 
except in high angle-of-attack operations, rotorcraft often exhibit undesirable 
pitch-roll coupling due to aircraft angular motion. For example, in response 
to a roll rate to the right, the TPP tilts to the left to provide the desired 
roll damping of the aircraft; however, as is evident from equation (13b), the 
TPP plane can also accompany a tilt in the for-aft direction, which produces 
an undesirable pitching moment. Experiment (ref, 5) has indicated that this 
kind of pitch-roll coupling can adversely affect helicopter flying qualities, 
especially in demanding tasks such as in NOE or IFR operations. It is desir- 
able, therefore, that the nature of the coupling in flapping due to aircraft 
angular rate be examined and a means of reducing the coupling be studied. 

First, we shall examine the effect of the primary design parameters on 
the coupling. Consider the hover case. From equations (13b) and (13c), the 
steady-state response of the longitudinal and lateral flapping to the aircraft 
angular rate in pitch and roll is given by 
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lA^ 

(p2 - 1)2 + (y 2/4)[(1/4) - (2/3)e + (e2/2)]2 



(17) 

In particular, when e = 0 and P = 1, equation (17) reduces to the following 
well known result: 



It is clear from equation (17) that either a pitch rate or a roll rate input 
will produce steady-state responses in both longitudinal and lateral flapping 
unless the following condition is satisfied: 


(y^/ 4)[(1/4) - (e/3)][(l/4) - (2/3)e + (e^/2)] 
2[1 + (eMg/Ig)] 


(19) 


where 



When the decoupling condition equation (19) is met, then 
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and 


Ba-i 


s. s. 


3b 1 
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9P, 
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- l/n[(Y/4) + (Y/8)(P2 - 1)] 


( 20 ) 


(p2 - 1)2 + (y 2/4)[(1/4) - (2/3)e + (e2/2)]2 

The physical meaning of the decoupling condition is given in the following: 

Consider the steady-state response of the longitudinal flapping ai due 
to aircraft roll rate p^. From equation (17), 


ai 


Pw/r2 


s.s. (p2 _ 1)2 + (^2/4) [( 1 / 4 ) _ (2/3)e +(e2/2)]2 

eM 




- 2(P^ -1)1+ 



( 21 ) 


Note that bi^ ^ due to has the identical expression. From the deriva- 

tion of the flapping equation, it can be seen that the first term in brackets 
in equation (21) is due to aerodynamic source while the second term is due to 
inertia. Thus, for P = 1 (the flapping frequency being 1 per rev), the iner- 
tia term vanishes, leaving only the aerodynamic contribution; in particular, 
for e = 0 and P = 1, ax = p^/iT^, as shown earlier in equation (18). How- 
ever, the flapping frequency may be changed appropriately by varying the hinge 
restraint, pitch-flap coupling, and hinge offset to achieve the level of iner- 
tia effect that offsets that due to aerodynamics. The condition for the iner- 
tia effect to cancel perfectly with the aerodynamic effect is essentially what 
equation (19) physically means. 

It may be instructive to examine in some detail some special cases of the 
decoupling condition, equation (19), and the associated flapping responses, 
equation (20). 

1, Rotor system with no flapping hinge offset— With e = 0, the decou- 
pling condition reduces to 


or 


P^ - 1 



128 


(22a) 


128 


(22b) 
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and the associated steady-state flapping becomes 


3ai 9bi 

s . s. s. s . 


9q 


w 


9p, 


w 


16 


( 23 ) 


It is interesting to note that the condition for decoupling the flapping 
responses due to aircraft angular rates p and q for the center-hinged rotor 
cannot be met without using either the pitch-flap coupling or hinge restraint, 
or both, as is evident from equation (18). Furthermore, when the decoupling 
condition is met, the flapping responses are identical to those for P = 1, 
irrespective of the values for K, and K . 

P 


Figure 15 shows the pitch-flap coupling required to decouple the steady- 
state flapping response due to aircraft angular rate, p and q, for a center- 
hinged rotor system. The required pitch-flap coupling is shown as a function 
of the Lock number and the hinge restraint. As can be seen from this figure, 
the value of the pitch-flap coupling in the sense of reducing the pitch for 
positive flapping (upward) increases as y increases and decreases somewhat 
as K increases. 

p 

2, Rotor system without flapping hinge restraint — With = 0, equa- 


tion (23) becomes 







(y^/ 4)[(1/4) - (c/3)] [(1/4) - (2/3 )e + (e^/2)] 
2[1 + (eMg/Ig)] 


(24) 


Consider that the rectangular blade has a uniform mass distribution. Then 
eM 3 /l 3 = 1.5 e. The pitch-flap coupling required to decouple the steady-state 
flapping response to the aircraft angular rate is obtained from equation (24) 
as follows: 


_ Y 1 - (8/3)e + 2e^ 12 e 

^16 1 + (3/2)e Y 1 - (4/3)e 


(25) 


Figure 16 shows the pitch-flap coupling required to decouple the steady-state 
flapping due to aircraft angular rate, p and q, for Kg = 0. The required 
pitch-flap coupling is shown as a function of the Lock number and the hinge 
offset. Note that for e = 0, it is identical with the curve in figure 15 for 
K 3 /lgfi^ = 0. However, as hinge offset increases, the required pitch-flap cou- 
pling differs substantially from that for hinge restraint with equivalent 
hinge offset (e.g., K 3 /lgfi^ = 0.225 has an equivalent hinge offset of 
e = 0.15, as far as flapping frequency is concerned). 

Figure 17 depicts the amount of hinge offset required to decouple the 
flapping response for a set of pitch-flap couplings. Note that for Kx = 0, 
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the amount of hinge offset varies almost linearly with Lock number from 
e = 0.019 for Y = 2 to e = 0.145 for y = 8. 

Aside from choosing a set of appropriate design parameters to decouple 
the steady-state flapping response due to aircraft angular rate, as discussed 
above, a feedback control law may also be used to achieve the same objective. 
Consider a control law of the form 



(26) 


which feeds back aircraft angular rate, p and q, to both longitudinal and lat- 
eral cyclic. The feedback gain matrix is to be chosen to achieve the 

desired decoupling characteristics: 



s. s. 



(27) 


In equation (27), the parameter ks is a desired sensitivity. Using equa- 
tions (13b) and (13c), it can be shown that the decoupling control law that 
achieves the decoupling characteristics of equation (27) is given by 


^ ^1 

*^f y[( 1/8) - (e/6)] 




+ k3(P^ - 1) 


( 28 ) 


2 

Q 




£ 

3 


+ 



Note that the sensitivity parameter in equation (27) is closely 

related to the contribution of the main rotor system to the vehicle’s damping 
in pitch and roll. Therefore, equation (28) provides an insight into the 
relationship among the decoupling control gains needed, the vehicle’s damping 
in pitch and roll, and the rotor system parameters. 
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CONCLUDING REMARKS 


The flapping equation of motion for a main rotor system that consists of 
four primary design parameters — flapping hinge restraint, flapping hinge 
offset, blade Lock number, and pitch-flap coupling — has been developed in this 
paper. Both teetering and nonteetering cases were considered in the develop- 
ment and the equation was presented in both rotating and nonrotating coordinate 
systems. 

Numerical examination of the tip-path plane dynamics indicated that the 
frequency of the regressing flapping mode was within the pilot’s effective 
frequency range (less than 2 Hz) when the Lock number was small (around 6 or 
less), especially when the hinge offset was large (s > 0.07). It was recom- 
mended, therefore, that the effects of flapping dynamics should be considered 
in assessing the flight dynamics and handling qualities of hingeless rotor 
helicopters or helicopters with heavy rotor blades. 

Study of the steady-state flapping response due to aircraft angular rate 
in pitch and roll at hover showed that there exists a condition for achieving 
a perfect decoupling, that is, roll rate results only in lateral flapping and 
pitch rate only in longitudinal flapping. The condition should prove useful 
in selecting appropriate values of the four design parameters to minimize 
undesirable pitch-roll coupling of helicopter flight dynamics, especially in 
demanding tasks such as nap-of- the-earth flight. A control law that feeds 
back aircraft angular rate in pitch and roll to both longitudinal and lateral 
cyclic can also be used to decouple the steady-state flapping response. It 
provides an insight into the relationship among the decoupling control gains 
needed, the vehicle’s damping in pitch and roll, and the rotor system 
parameters . 

The cyclic control phasing or control mixing necessary for decoupling the 
steady-state response of the longitudinal and lateral flapping has been devel- 
oped for hover and for forward flight. Because of its simplicity, the decou- 
pling control law for hovering flight has been used previously in a ground 
simulation. It has provided an accuracy adequate for the entire low speed 
flight regime. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. 94035, May 11, 1979 
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APPENDIX A 


DERIVATION OF FLAPPING EQUATION 

Using the assumptions described in "Flapping Equation of Motion" (p. 2 
of this report), the flapping equation was derived by summing the moments at 
the flapping hinge. 

"a + »CF + “l + “cor. + “r + “bA + “bL + “w ■ 0 

where 

moment due to aerodynamic force acting on the blade 
moment due to the centrifugal force, cos 3 + eM^]sin 3 

moment due to blade inertia, = -Io3 

■L P 

moment due to Coriolis acceleration, 

Mcor. = 2[I^ + eMg](pS^ cos ip' - sin ip/ ) 

restraint moment, = “Kg3 

moment due to body angular acceleration, “ ^3^^ 

moment due to body normal acceleration, = M^(w - uq + pv) 

weight moment of the blade about the flapping hinge 





The moment due to the aerodynamic force acting on the blade, will now be 

derived in the following: 
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At an azimuth station is given by 




J rR-e 

D 


dF r* 


where 


dF^ = (QR)2ac[Ux^0 + U^Upldr’ 


and 


U = 

T — = e(l ~ cos 3) + y sin + x cos 3 


U _ r/ 

Up 4 = X cos e - y sin B cos ij; - S(x - e) + X cos 3^ + -^ sin 


- sin B._ - cos 


w Q 


sJcosK.] 


y ^ 


A V cos g 


QR 


. V sin a - v^ 

\ ±_ 

QR 

A e + r’ 


R ’ ^=t’ 


0 = 00 - A cos 1/^-3 sin Jp + x0^ - K, 3 
ic ic t 1 

For small 3, 

U-p = X + p sin ij; 

Up = A - yg cos t/j - 6(x - e) + x|j^— cos 3^ + ^ sin B^^sin \p 


■ (n ^ ^w) 


and may be found to be 


sin 
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^ (1 _ e**) + i (1 - e^) (2y sin .4- - e) + 4 (1 - e^) (y^ sin^ 4' 

(p/2)ac(fiR)"R^ 3 ^ 

- 2ye sin 4) - (1 - e)y^e sin^ 4 jJ (6o - cos ^ sin 4) 

+ (1 - e®) + -^ (1 - e'*)(2y sin 4 - e) 

+ -i (1 - e^)(y^ sin^ ip - 2ye sin 4<) - y (1 - e^)y^e sin^ 4<J 9^ 
+ jj (1 - e^) + j (1 - e^)(y sin 4' - e) - (1 - e)ey sin A 


- ( [i (1 - + i (1 - 

+ -| (1 - e^)(u^ sin^ 


e^) (2u sin ip - e) 

- 2pe sin \p) - (1 - e)u^e sin 
e^)(]i sin i|j - e) 


2 


- (1 - e)ue sin ip 


j y cos ^ 




- (1 - e**) + -| (1 - e^)(y sin ip - 2e) 

- (1 - e^)(2ye sin 4 - e^) + (1 - e)ye^ sin i/; j B 

+ (1 - £**) + "2 (1 ~ £^)(y s in ip — e) 

- ^ (1 - c"^)sin j(§ cos 6^ + ^ sin B^)sin 4 

- sin 3 - 4 cos B ^ cos 4^1 ^^2) 

\Q w Q w/ 

Neglecting the terms containing e^, e**, . . the following result is 

obtained : 
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(p/2)ac(nR) 


" ■ f + sin - e + -|u sin - eu sin + -| ye^ sin 

^ (So - cos tp - Bjj, sin ij;) + [^ ~ |- + U sin ^ ~ "I s 
+ j M sin - I ue sin i e + y sin ip 

(l ■ ^ ^)] ^ ~ {[i ~ I sin (I - e + ^ sin 

ey sin p + — y sin ~ ^ ^ ^(~t ~ ^ 


y cos ip 


€ -h C' 


+ y sin — 1 


— sin B, - -^ 


w n 


^ w ) 


The flapping equation is then obtained by substituting all the moments into 
equation (Al). The result becomes 


1 H-* 


^ [i ■ I + T + ^ + T 

. Ifl e /2 y . p2 . 

2 1 1^4 3"^^\3 c + sin y sin ijy'j sin 

= 2 ^1 + J (pQ cos i|j' - qQ sin i{/ ' ) + (p sin ip ' + q cos i|;') 


).] 


y , e ^ 

^ "2 ^ ~ sin ip + y sin ijy' 


y cos Ms 


^ sin 4.j 


= 2(1 + 


4 . p r 
+ — (w - uq + 


pv) -^ + ^/ ri _£ 

^3 2 \U 3 


+ y sin ^ 


e + y y sin 


- zp sin P + ^ yc2 sin ^)] (0„ - cos 4> ~ sin ^) + [^1 - | + y 

X (I - I e + ^ y sin i|- - I ye sin ((;) sin 4/j 6^ + [^^ - i e + y 

"" G ■ ^ I '^J ^ 1 ■'■ Ki “ f) 

^ (Continued) 
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- sin - qfi cos B^) cos 'I'J + ^ p(f ~ c) cos B^ + sin B„) 

X (1 - cos li>) - (p{^ sin B^^ - qf5 cos Bw)sin 
where the Lock number y is defined by 


(A4) 

(Concluded) 


Y ^ 


A P^cR^ 
‘B 


Equation (A4) may be rewritten in a nondimensional form by defining 


i A 1 B = ^ 

^ ^ di^ 


1 d"B 

di|.'^ 


A B = 


The result becomes 


l+il4--le + ^ + 


+ 




[i- 

[4 3 \3 


(I - s sin (.jsi- j^l + 

U sin i|,'^u sin 'JjJ 


E + -^ sin ij,' - ep sin + 


[i- 

c- 

— "i 1 

/I , c \ ■ 

It: - e+^I y sin y 

y c 0 s i|i 1 

L3 

2 

\2 If 

) 




COS - t: sin ijj 


M. 




(w - uq + pv) 


M, 


w 




lofi" 2 L4 3 


sin u - f:u sin v 


-I- ~ ye^ sin 


in i^'^y sin iJcSo “ ^ic ” ®ic [^5 4 

+ 1 iT + 1 u sin li,. - X ys sin u sin tlj “ I ^ 

+ (i - E + I c') P sin ;] x| +7(^-3) j^(t + o 

- (I sii^ “w ~ f "] ^ 8 ■ ")[( 

X (1 - cos 2i) - sin cos sin 2cj 


COS 8 +— sin 8 

w ' 


.) 


(A5) 
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APPENDIX B 


TRANSFORMATION OF FLAPPING EQUATION FROM ROTATING TO NONROTATING 

COORDINATE SYSTEMS 


In rotating coordinate systems, the flapping equation is given by equa- 
tion (1) and is repeated here. 


’^i ■ 3 ^ ^ '^i\ 6i+ f) 

V (i - " + x) 2^i] + V (I - ") ^ 

+ Y (i - ^ + t) (^ ■ 2^i)j I B. 

" f (“ '^i if ’^i) + ^ "^i + ^ 


COS ij; , 
1 


cos i) , 


Ig^ f - uq + pv) - g] + 5 1^(1 - f) + (I - e)p sin . 

^ ~ ^ - cos 0 0 

■ 2 [(i - f) ■=»“ *1 + “(l‘ f) (i - E + 

- cos - I [(i - f)=“ *1 + “(j - f) (* - 2*i) 

+ X (i ' ^ + t) *i - Bic + 2 [(i ■ f) 

+ c(| - I c) sin ♦. + (i - l) (1 - cos 2t^)j 6^ + 1 [(i - |) 

+ S0 - c + ^) sin 1 . ,j j| + f (I - f)(^ sin ^ cos 1 J 

8 "(I ■ ') L’j? * X 2*l] 


(Bl) 
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where 


eM^ 


K, 


= 1 + 








i = 1, 2, 3, . . N 
= p cos 6^ + q sin 
= -p sin Bw ^ ®w 


To gain a better insight into the dynamics of the flapping motion equa 
tion (1) will be transformed into a nonrotating coordinate system using the 
multibladed coordinate transformation (ref. 9), that is. 


where 


nc 


ns 


k 


= Bo + + X) 

Pnc 

+ 

^ns 

sin nij^^ 

n=l 

i = 

1. 2, 

3, 

. . ., N 


1 

2 

(N - 

1) 

N odd 

1 

2 

(N - 

2) 

N even 


N 



1 

“ N 

E 

^i 



i=i 



1 

N 

N 



E 

B.(-l)^ 

(Bj = 0, if N odd) 


i=l 


= I E 

i=i 
2 ^ 

= N E '"'^1 

i=i 

. ♦ (1 - 1 ) 


(B2) 


Let Bo and B be denoted, respectively, by 
— R — 

^ i (Bi,B2 3n)^ 

_6 = (3o’^ic*Pis’^2c’^2s» • ’ •’ ^(1/2) (N-l)c’ ^(1/2) (N-l)s^ 

T 

= (6o,Bjc>^is ®(i/2)(N-1)c’ ^(1/2) (N-l)s’ ®d^ 


for N = odd 
for N = even 
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Then the transformation equation (B2) may be rewritten by 


where, for N = odd 
r 1 cos 1^1 sin 


cos Y (N - sin y (N - 


1 cos i '2 sin Tpz 


cos y (N - Sin 4 (N - l)i ['2 


1 cos \p . sin ^1' . 

1 1 


cos ^ (N - sin ^ (N - l)ii 

^ 1 Z i 


1 cos sin ip^ 


and, for N = 


cos j (N - l)ip sin ^ (N - l)ip 


cos 

sin 

cos (N - 

2)^1 

sin 

1 (N - 

2)^1 

(-1) 

cos 

sin ip 2 

cos i (N - 

2)<p2 

sin 

j (N - 

2)]p2 

(-1) 


1 cos \P‘^ sin 


cos Y (N - 2)t|^i sin y (N - 2)ip^ (-1)^ 


1 cos sin 


1 cos sin 

. N N 


cos - (N - sin j (N - (-1)^"' 


cos j (N - 2)^^ sin I (N - 2)^^ (-1)^ 


"^'i “ + IT ’ i = 1,2,3, . . ., N 
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The flapping equation (Al) may now be rewritten in a matrix-vector format 
as follows: 

(B4) 

where 



f^ a right-hand side of equation (B-1). (B4) 
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With respect to the nonrotating coordinate system the flapping equation 

(B4) may be represented by (B5) using the transformation (B3) . 

B. + D_6 + KB = f^ (B5) 


where 


i ■ 

D = T“’^(2T + AT) 

K = T“^(T + AT + BT) 


To expedite the calculation of the Inverse of T, a proper scaling of the 
variables may be made so that 

I = TgB.j (B6) 


where 



0 


for N odd 


0 


T 


S 


r 1 


/N 




0 


0 



for N even 


The resultant transformation Tj = TTg is orthogonal. Since 

the desired simplification in computing T 


is achieved. 


Tables 1 and 2 show the results for N = 3 and N = 4, respectively. For 
a special case in which Ki = e = 0 and p=q=p=q=0, the equations in 
these tables reduce to those previously obtained in reference 8. 
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APPENDIX C 


FEEDFORWARD CONTROL LAW FOR DECOUPLING FLAPPING RESPONSE 

IN FORWARD FLIGHT 

From table 3, it is evident that in forward flight the coning is coupled 
to the cyclic flapping. To achieve a perfect decoupling, a complete control 
mixing will be required 

Consider a feedforward control law of the form 


"00 " 


r<5 1 

c 

A 

= K., 

6^ 

is 

D 

A 

JlS- 




(Cl) 


where 


V A 


36 


36. 


300 

300 

900 

36 

c 

36^ 

36l 

'^s 



9«c 

36a 

36l 

3Bls 




36. 


is the "control mixing" matrix. The control law is to be selected so as to 
achieve a desired decoupling in steady-state response having the sensitivities 


1, A ^^0 
0 = 86 


ki 4 

^ 86 . 


and ko ^ 


A 


86t 


as follows: 


“ao " 


a. 


Is 


_^s. 

s . s . 


ko 

0 

0 


0 

0 

k„ 


0 

ki 

0 


>-^L' 


(C2) 
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Denote 


a, = 


Sp ■ <«c’‘a’V" 



Let the control effective matrix E be 



and 

is, 


and ^2 account for the effect of inflow due to collective input. 



that 


The steady-state response of the coning and cyclic flapping to control input 
is given by 
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a = K“^f 

— s . s . — 

= K'^Eu 
= 

Imposing the desired decoupling characteristics as described in equation (C2) , 
that is, g = SUp, the control law may be found to be 

Kp = E"^KS 

It is readily verifiable that 

"“5 

, 1 

^ 0 

Cl {a^a^ - a^a^) 

L-“4 

and the desired decoupling control law is obtained as described in equa-- 
tion (C4). 
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TABLE 1.- FLAPPING EQUATION IN NONROTATING COORDINATE SYSTEM, N = 3 


I + D3i + K3B = f 3 ; i fi (Bo.8ic.6ls) 


0 
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i rp(| - + c^) 

2 - ^ - s: + e^^cos 3ij( 


Ks 3 4 ) ^ (3 
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¥ (i - ■-) 
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TABLE 1.- Concluded. 
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TABLE 2.- FLAPPING EQUATION IN NONROTATING COORDINATE SYSTEM, N = 4 


B + D 4 B + K 4 I = I 4 ; 3 A 


+ ^) 
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TABLE 2.- Concluded 



37 


TABLE 3.- TIP-PATH PLANE EQUATION 


a + Da + Ka = f 
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TABLE A.- COMPARISON OF CHARACTERISTIC ROOTS OF THE REDUCED AND UNREDUCED 

TPP MODELS AT y = 0.3 
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Figure 1.- Main rotor primary design parameters. 




Figure 2.- Coordinate system for the main rotor. 
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Figure 4.- Effect of hinge restraint and Lock number on TPP modes at hover. 
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(jO/Q> 



Figure 5.- Effect of pitch-flap coupling and lock number on TPP modes 

at hover. 
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Figure 7.- Effect of advance ratio on the TPP dynamics, e = 0.15, K 












TIME, sec 


(d) Flapping. 
Figure 13.- Concluded 







TIME, sec 


(a) Coning, 



0 .2 .4 .6 .8 1.0 

TIME, sec 


(b) Longitudinal flapping. 

Figure 14.- Comparison of responses of the reduced and unreduced order TPP 
dynamics to unit step input in longitudinal cyclic (K^ = Kg = 0, e = 0. 
y = 0. 3, Y = 8 






-25 


-50 

( 

LOCK NUMBER, 7 

Figure 15.*- Pitch-flap coupling required for decoupling flapping due to 
aircraft angular rate for center-hinged rotor at hover. 



Figure 16.- Pitch-flap coupling required for decoupling flapping due to 

aircraft angular rate at hover. 






2. Government Accession No. 


3. Recipient's Catalog No. 


1. Report No. 

NASA TP- 1431 

4. Title and Subtitle 

EFFECTS OF PRIMARY ROTOR PARAMETERS ON 
FLAPPING DYNAMICS 

7. Author{s) 

Robert T. N. Chen 

9. Performing Organization Name and Address 

NASA Ames Research Center 
Moffett Field, Calif. 94035 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

15. Supplementary Notes 


16. Abstract 

An analysis has been made to study the effects of flapping dynamics of 
four main rotor design features that influence the agility, stability, and 
operational safety of helicopters. The parameters include flapping hinge 
offset, flapping hinge restraint, pitch-flap coupling, and blade Lock 
number. First, the flapping equations of motion are derived that explicitly 
contain the design parameters. The dynamic equations are then developed for 
the tip-path plane, and the influence of Individtial and combined variations 
in the design parameters determined. 

The analysis includes a study of the steady-state flapping response 
with respect to control input and aircraft angular rate which leads to a 
feedforward control law for control decoupling through cross feed, and a 
feedback control law to decouple the steady-state flapping response. The 
condition for achieving perfect decoupling of the flapping response due to 
aircraft pitch and roll rates without using feedback control is also found 
for the hover case. 

The analysis also indicates that the frequency of the regressing flap- 
ping mode of the rotor system can, for some designs, become low enough to 
require consideration in the assessment of handling characteristics (less 
than 2 Hz). This occurs for rotors with a large effective hinge offset and 
with heavy rotor blades. 
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